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Abstract 

We consider periodic traveling wave solutions to the focusing nonlinear Schrodinger 
equation (NLS) that have been shown to persist when the NLS is perturbed to the 
complex Ginz burg-Landau equation (CGL). In particular, we show that these periodic 
traveling waves are spectrally stable solutions of NLS with respect to periodic pertur- 
bations. Furthermore, we use an argument based on the Fredholm alternative to find 
an instability criterion for the persisting solutions to CGL. 

1 Introduction 

The complex Ginzburg-Landau equation (CGL) is used in several contexts such as chemistry, 
fluid dynamics, and optics [H HTJ HE]. Furthermore, CGL was shown to be the generic 
equation modeling the slowly varying amplitudes of post-critical Rayleigh-Benard convection 
[22] . CGL is not in general an integrable equation. However, a successful strategy when it 
comes to the study of existence and stability of solutions is to consider CGL as a perturbation 
of the nonlinear Shrodinger equation (NLS). One can then take advantage of the fact that 
NLS is completely integrable. The drawback of that method is that the coefficients of the 
dissipative terms in CGL must remain small. 

In the case of localized traveling wave solutions of CGL such as pulses and holes, this 
strategy has been implemented in several instances (see for example [151 HSl [IS])- the 
case of periodic traveling waves, it has recently been shown that classes of solutions of NLS 
persist when considering CGL perturbations [HI |9l HO] . The main purpose of this paper is to 
study the stability of one of the classes of such solutions, namely the cnoidal wave solutions. 
We will restrict ourselves to the study of stability with respect to perturbations that have 
the same period as the waves which amounts to restricting the spatial domain of the solution 
to an interval of length equal to the period. 

Our results are twofold. First, we describe the spectrum of the linearization of the NLS 
about these cnoidal waves defined by ^ below. In particular, we show that, when considering 

* Supported by NSF grant DMS-0608587 
t Supported by NSF grant DMS-0509622 



1 



perturbations that have the same periodicity as the wave, the cnoidal solutions are spectrally 
stable, i.e., the point spectrum of the linear operator arising from the linearization do not 
have eigenvalues with strictly positive real part. In fact, we show that the location of the 
point spectrum is restricted to the imaginary axis. In addition, we obtain the continuous 
spectrum and show that it intersects the right side of the complex plane, showing that cnoidal 
wave solutions ^ are unstable when considered on the whole line. On the one hand, this is 
consistent with the fact that complex double points in the spectrum of the Lax pair usually 
give rise to instabilities fHX [T21 120] ■ On the other hand, it is also consistent with recent work 
^ in which ([s]) is numerically shown to be unstable. 

Our results concerning the spectrum of linearized NLS are summarized in Theorems 
[l] and [2] and illustrated in Figure |5| The main tool we use is the fact that solutions to 
linearized NLS can be constructed using squared AKNS eigenfunctions [121 EI] ) ^oi which 
explicit formulas are available [31 IS]- The main difficulty lies in the fact that one has to 
show that all the solutions of the eigenvalue problem are obtained that way; in fact, this 
is not true for the zero eigenvalue. By using an alternate route to construct solutions, and 
employing the technique of the Evans function (as defined in [I3|), we show that the zero 
eigenvalue has geometric multiplicity. 

Secondly, we use the knowledge of the spectrum of the linearized NLS to carry out a 
perturbative study of the spectrum of linearized CGL. More precisely, we linearize about 
the persisting periodic solution of CGL, and we use the Evans function technique and the 
Fredholm alternative to identify criteria for linear instability. Under these conditions, eigen- 
values of the linearized CGL bifurcate from the origin to the right side of the complex plane 
as NLS is perturbed to CGL. These results are summarized in Theorem |3] and illustrated in 
Figures [T] and [2j 

The article is arranged as follows. In ^ we define the concept of Evans function for 
the linearized NLS about a periodic traveling wave. We then determine the spectrum of the 
linearized NLS about the cnoidal waves. In ^ we describe the persisting solution of CGL 
and the corresponding linearization. In ^ we use the concept of Evans function and the 
Fredholm alternative to study eigenvalues emerging from the origin when NLS is perturbed. 
The results of more detailed calculations are relegated to Appendix |X] and the proof of 
Theorem |2] is given in Appendix [Bj 
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2 Evans function for the linearized NLS about a peri- 
odic traveling wave 

The NLS equation is written in system form as 

igii + qixx + 2gig2 = 
-mt + q2xx + 2gig2 = 0. 
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e U2, for a real constant a, and obtain 



(2) 



The reality condition q2 = qi gives the focusing NLS for q = qi- 

We make a change of variables qi = e~'°*Mi, q- 
the following system: 

iuit + aui + Uixx + '2,u\u2 = 

— iu2t — au2 + U2xx + 2uiul = 

This system has several T-periodic traveling wave solutions of the form ui = U2 = Uq{x — ct) 
where Uq is expressible in terms of elliptic functions. However, these only persist as T- 
periodic solutions of the complex Ginzburg-Landau (CGL) equation when c = [U [9l [TO] . 
Thus, we limit our attention here to the real solution of (|2]) given by 

t) = U2(x, t) = f/o(x), f/o(x) = 5A; cn(5x; fc), (3) 

with a = (5^(1 — 2k'^) where k is the elliptic modulus, < /c < 1. The period of this solution 
is T = 4K/6, where K = K{k) is the complete elliptic integral of the first kind. 

We then linearize the system ^ about the solution ([3]). Substituting ui = Uq + yi, 
U2 = Uq + V2 in ^ and discarding higher-order terms in the v's, we get 



ivij + avi + vi^^ + 4f/oVi 

-iV2t + aV2 + V2x 

We substitute into Q the ansatz 

vi = e^^Wi{x) 



2U^V2 

+ AU^Y2 + 2U^Yi 




0. 



V2 = e^%2(a;) 



(4) 



(5) 



where wi and W2 are assumed to be independent of time t, obtaining the following ODE 
system for wi, W2 

iiwi + awi + wixx + 4:UqWi + 2UqW2 = 0, 
-i£w2 + aw2 + W2xx + ^Ulw2 + 2Ulwi = 0, 
which can be written as the eigenvalue problem 



CoLJ 



iWi,W2) 



where 



Co = J 



dx" 



, , 4t/o^ + « 



2Ui 



J 







(6) 



(7) 



'0 ^"^0 

Since we are interested in periodic perturbations, we define the point spectrum of the eigen- 
value problem (j6| to consist of values of £ for which lo is an element of the space (L2(T))^ of 
two-dimensional vectors with components in the set of T-periodic square integrable functions. 

In order to define the Evans function, we rewrite the ODE system ^ as the following 
first-order linear system with T-periodic coefficients: 



dx 



w = A{i,x)w, w = {wi, W2, w[, W2] 



where A{i, x) is the 4x4 matrix given by 



A{i,x) 



( 



-(a + i£ + 4f/2) 
. -2^0^ 







-2Ul 



1 







(9) 



-{a-\i + MJl) 0/ 



Let be a fundamental matrix for this system, with $(0;^) = /. Then the Evans 

function [13] is defined to be 

M{k, i) = det (<I>(T; £) - e'^'l) . (10) 

If M(K,i) = for some k G M then the system (|8]), or equivalently (|6|, has a bounded 
solution. The zero set of M(k,£) thus corresponds to the continuous spectrum of Co- How- 
ever, we restrict the eigenvalue problem ^ to T-periodic square integrable functions which 
corresponds to setting k = 0. We thus define 

Mo{£) = M{0,i). (11) 

The zero set of Mo{i) corresponds to the point spectrum of the eigenvalue problem ^ and 
the multiplicity of the zero of Mq corresponds to the algebraic multiplicity of each eigenvalue 



For all but a few exceptional values of £, the Baker eigenfunctions associated to periodic 
traveling wave solutions of NLS may be used to construct a fundamental solution matrix 
for the system dsl). (Among the exceptions is the important case i = 0, which is treated in 



^4.1 ) Then the Floquet discriminant of the AKNS system may be used to determine the 
zeros of the Evans function. The following theorem, which establishes the spectral stability 
of the NLS solution ^ with respect to T-periodic perturbations, is proven in Appendix |B] 



Theorem 1. The zero set of the Evans function ( 11) which corresponds to the point spectrum 
of Cq defined in ^ consists of a countably infinite number of points on the imaginary axis 
in the i-plane. With the possible exception of four values, the locations of these points is 
determined by the formula i = 4i/i where 



/i2 = (A2-A?)(A2-Ai'), (12) 

and X is a value for which the Floquet discriminant A(A) of the AKNS system equals ±2 or 
0. 



Note that the branch point Ai is related to the solution (|3]) by Ai/|Ai| = —k' + ik and 
= 2 I All, where k' = \^T^^W is the complementary modulus and the Floquet discriminant 



is defined in (73). The exceptional values are those for which (12), as a polynomial equation 
for A, has multiple roots. Note that these all correspond to imaginary values of i. 

Figure |5] shows the location of the spectrum in the £-plane for 6 = 2 and three different 
values of k. In addition, the continuous spectrum is also obtained in Appendix B} It consists 
of the imaginary axis and a figure eight centered at the origin (see Figure [5). Thus, this 
shows that solution ^ is unstable with respect to bounded perturbations and is consistent 
with the fact that complex double points in the spectrum of the Lax pair usually give rise 
to instabilities [TT|[T2ll20]. 

3 The perturbed setting 

3.1 Periodic traveling wave solutions to CGL 

The CGL perturbation of the NLS system has the following form: 

kit + qixx + 2gig2 = e[irqi + iqi^^ - 2isqlq2\, ^^^^ 
-iq2t + q2xx + 2qiql = e[-irg2 - iq2xx + ^isqiqj]. 
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Here, e > 0, and r, s are real parameters. Again, g2 = ^ is the reality condition for the 
focusing case. 

We again make the change of variables qi = e~'"*Mi, q2 = e^°*M2, and obtain the following 
system: 

iuu + aui + uixx + 2u{u2 = e[irui + iui^^ - 2isulu2\ 

.2 



l\u2 = e[irui + iui^^ — ^isult 
-m2t - au2 + U2xx + 2uiul = e[-irM2 - m2xx + 2isuiU2] 
u. 



(14) 



with the reality condition U2 = ui 

The solution to NLS given in ([s]) persists as a T-periodic solution of the CGL (14) if r 
and s satisfy a certain algebraic condition. More precisely, it was shown in [HI El HO] that for 
small e > 0, the CGL equation (14) has a solution which is T-periodic in x of the form 

u{x,t) = U{x) with U{x) = Uo{x) + eUi{x) + 0{e^), (15) 



(16) 



where Uq is the solution of NLS given in (|3]), provided that r and s satisfy 

and r > 



3{E-Kk'^) 



where E = E{k) is the complete elliptic integral of the second kind. 

In order to investigate stability properties of the CGL solution U, we will need the 
expression giving Ui. To find it, we substitute ui = U and U2 = U into (14). At first order 

in e, one finds that the vector U = (^Ui, Ui, U[, Ui j satisfies a non-homogeneous linear 
dynamical system. This system is given by 



dx 



U = ^(0,x)U + N, 



(17) 



where .4(0, x) is obtained by setting £ = in the matrix given in ([9]), and the non-homogenous 
term N is given by 







N 



ir[/o + if/o -2is[/o^ 
-irUo-iU;; + 2isUl 



The homogeneous system corresponding to (17) is given by 

d 



dx 



U = ^(0,x)U. 



;i8) 



It turns out that (18) can be solved explicitly. To do so, one uses two symmetries of CGL 
(14): translation of the variable x and phase change —>■ g^e*^ for any real 6. The presence 
of these two symmetries implies the following T-periodic solutions to (18): 



(19) 









'Uo' 




K 


yh _ 


-Uo 
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One can then use these two solutions to perform reduction of order on the homogeneous sys- 
tem. The two-dimensional homogeneous linear dynamical system for two unknowns obtained 
that way turns out to be diagonal and thus can be solved by quadrature. The expressions for 
the other two solutions Vg and are rather complicated. They are given in Appendix|A| No 



linear combination of V3 and V4 is periodic, implying that (|18l) has only a two-dimensional 



space of T-periodic solutions generated by and Vj. Once four linearly independent so- 



lutions to (18) are known, one can then use the variation of parameters method to find a 



particular solution U to (17). Requiring this solution to be T-periodic leads to the algebraic 



relations (16) between the constants k, r, 6, and s. The first component of U is then given 

(20) 



Ui = iSk{s + l)cn(fe) 



+ 



K (2 dn{6 x) sn{6 x) Z{5 x) - Z\6 x) cn{6 x) + cn^{6 x) k 



2\-\ 



3 cn(5x)(E + (P - i)K) 

where Z is the Jacobi zeta function and 64 is a Jacobi theta function. (For definitions and 
conventions for these and other relatives of the elliptic functions, we refer the reader to [1].) 

Note that any linear combination of or can be added to U without losing the peri- 
odicity property. This latter fact can be explained by the presence of the two symmetries of 
CGL mentioned above: adding a multiple of the top component of V^* corresponds to adding 
translation in x to the CGL perturbation, while adding a multiple of the top component of 
V2 corresponds to adding a phase change. 

Now that we know more about the dependence of U (x) in terms of e, we can study its 
spectral stability as a solution of CGL 



3.2 Linearization of CGL 



We linearize the system of equations (14) about the stationary T-periodic solution ui{x, t) 



U{x), U2{x,t) = U{x), where U is given in (15). We do this in a similar way as it is done for 



NLS in §2l Setting ui = U + vi, U2 = U -\- V2 and discarding all higher-order terms in the 
v,'s, we obtain 



(21) 



(22) 



ivu + avi + vi^.^ 4|f/pvi + 2U'^V2 = e[irvi + ivi^.^. - 4is|?7pvi - 2isf/^V2], 
-iv2t + av2 + via-a; + 4:\U\\2 + 2f/Vi = e[-irv2 - iv2xx + 4is|[/|\2 + 2is[/^Vi]. 

Next, by analogy with [|2]we substitute into ( [2I] ) the ansatz 

vi{x,t) = 6^*1^1 (a;), V2(x,t) = e^^W2{x). 

This results in the following ODE system for Wi, W2 

I I O Or I I O On 

iiwi + awi + wixx + 4|f/| wi + 2U W2 = e[irwi + iwixx — 4is|t/| wi — 2isU W2\ 
-iiw2 + aw2 + wixx + 4|[/pW2 + 2U'^wi = e[-irw2 - iW2xx + 4is|f/|^w;2 + 2isU'^wi], 

which can be rewritten as the eigenvalue problem 

LijO = iijJ^ a; := (wi, ^2)^ . (23) 
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If we expand L = Co+^'^i+^^'^2+0{e^), then Cq is the operator arising from the hnearization 
of NLS given in ([T]) and Ci is given by 

r-T^,f r~AsU^ 2Uoi2iU^-sUo)\ . _ A 0\ 
' 'dx2 \2Uo (2if/i - sUo) r - 4sU^ J ' VO V ' 



Since we are interested in periodic perturbations, we restrict the eigenvalue problem (23) to 
the space (L2(T))^ of T-periodic square integrable functions. 

4 Eigenvalues near the origin 

In this section, we study stability properties of U (x) as a solution to CGL by studying the 
spectrum of the operator C. We use a perturbative argument to study the eigenvalues of £ 
near i = 0. In order to do this, we first need to characterize the eigenvalue £ = for the 
problem ^ arising from the linearization of NLS about the T-periodic solution ([s]). 

4.1 The case e = 

In this case, we show that the algebraic multiplicity of the zero eigenvalue is four by showing 



that the lowest order nonzero derivative of the Evans function (11) at £ = is the fourth 
derivative. More precisely, we have the following theorem 



Theorem 2. The Evans function Mo{£) defined in (11) is such that 



Mo(0) = M^(0) = K(o) = K'(o), 

where the primes denote derivatives with respect to i. Furthermore, the fourth derivative, 
which is given by 

is nonzero for all k G (0, 1). Thus, the eigenvalue i = is of algebraic multiplicity exactly 
four. 



Proof. Let Wi{x;i), i = 1 ... 4, be four linearly independent solutions of (|8]) and let W{x; i) 
be a fundamental matrix solution of ([s]) with columns Wj. In order to compute the derivative 
of the Evans function at £ = 0, we expand each of the solutions as 

Wi(x; £) = w°(x) + iwlix) + fw^{x) + 0{e). (25) 

The w°(x) are solutions of (pi) for i = 0, so they are linear combinations of the , i = 1 ... 4 



given in ( 19 ) and in Appendix |A| We make the choice 

w°(x)=Vf(a;). 

Recursively, the higher order terms for j > 1 satisfy the non-homogeneous system 

^w^=^(0,x)w^ + Mr\ (26) 
7 



where 




and (w^)^ denotes the kth component of the vector w^. Since we have a fundamental set of 
solutions for the homogeneous system, we can solve (26) by variation of parameters method. 
For the higher order terms we choose 



Wo{x) / W,-\z)M{{z)dz, 



1...4, j>l, 



where Wo{x) = W{x; 0) is the fundamental matrix solution for the system ([8]) at 
columns V^. The formulas for w} and can be written in a compact way as 



with 



/ -xUo 
xUq 
- {xUo)' 



Wo 



( (xt/o, 
X Uo)" 



(27) 



where Uo is the solution of NLS given in ([s]). The expressions for w|, W4, wf, and W2 are 
given in Appendix |Aj 

In order to compute the Evans function (11), we need solutions i = 1 . . .4 of 

(p| with initial condition 0^(0, = Sj, where ej are the standard basis vectors of M^. The 



i) are thus the columns of the matrix $(a;; i) used in the definition of the Evans function 

dth those 

(W3 + W2) 



(10). Matching the latter initial conditions with those of the Wj, one finds that 

1 



01 

02 
03 
04 



2Sk 
1 

25k 
1 

26H: 
-1 

Wk 



(Wg - W2) , 
(5^W4 — Wi) 
(5^W4 + Wi) 



The Evans function (11) can then be expressed as 

1 



Mo{i) 



Mo(£) 



where 



Mo{i) = det {W {T;i) - A) , A 






-S^k 
5^k 



5k 5k " 

-5k 5k 

5k 

5k 



and W{x;£) is the fundamental matrix solution of ([8]) with columns Wj, i 
that W{Q;e) = A. 



{2t 



(29) 



1...4. Note 
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Since Mo(£) differs from Mo(£) by a factor wliicli is constant in ^, tlie derivatives of Mq at 
i = are constant multiples of those of Mq. The derivatives of Mq at £ = can be computed 
using the expansion (25). Because Vf,? = 1,2 are T-periodic, the first two columns of 
W{T] t) — A are zero at £ = 0. Thus, the expansion of the determinant Mq in powers of £ 
has no terms of lower order than C.'^. Furthermore, the coefficient of ^^ is given by 



-<(0) = det (w}(T), w^(T), V3^T) - A,,\l{T) - A,) 



(30) 



where Ai denotes the iih column of A. 



The right-hand-side of (30) is easily shown to be zero, as follows. Using the expressions 
for w\ and V4 given in (27) and in Appendix [A| we find that the vector 



w 



9,1 



(31) 



is a T-periodic function of x. Taking the difference of the values of the vector w^^i at x = 
and X = T gives a linear combination of the first and fourth columns of the matrix in (30) 
which must be zero. Similarly, the vector 



w,,2 = W2 - ^i? {k"K + E{1 - 2^)) V3^ (32) 

is periodic and thus a linear combination of the second and third columns of the matrix in 
(30) is zero. 

These facts also imply that Mq"(0) = 0, as follows. The idea is to write Mo(£) as 

^ = det (w}(T) + 0(£), w^(T) + 0(£), V3^(T) -A, + 0(£), Vj(T) - A, + 0{l)) . 

(33) 

Since the term of the right-hand-side of (33) will be a sum of determinants of matrices 

ero 
0. 



each involving three of the four columns of the matrix in (30), it must be zero. Thus, we 
have that the Evans function for e 
Then, 



at least has a fourth order zero at 



Mo( 



+ 
+ 
+ 



det [ wl{T) + wl{T)£ + Oif), wl{T) + wl{T)£ + Oif), 

V'^{T) -A3 + wl{T)i + Oif), V^^{T) - + wl{T)i + Oif 
det{wl{T),wl{T),y^,iT) - AsMm - A,) 
det( w2(T),wi(T),wi(T),V4^(T) - A4) 
det( w}(T),w2(T),V3^(T) - A3,wl(T)) 
'deti wliT),wliT),wliT),wliT)) + Oif). 



(34) 



It is then straightforward to use (34), and the expressions for w,^, w^, and in (27) and 
in Appendix lAl to find 



M(S^)(0) = 1536 



(^2 + K{K - 2E)k 



(35) 



Then one obtains (24) using the relation (28). Using the inequahties K > E > k'K, which 



hold for all moduli k G (0, 1), we find that 

E"^ + K{K - 2E)k'^ > 2k'^K{K - E) > 0, 
and thus the fourth derivative of the Evans function is nonzero for all such k. 



□ 



We have established in ^3.1 that for for £ = (|8]) only has a two-dimensional space of 
T-periodic solutions generated by Vf' and Vg. The eigenvalue problem ^ thus admits a 
two-dimensional eigenvector space for i = 0. However, since the zero of the Evans function 
is of multiplicity four, the eigenvalue i = actually is of algebraic multiplicity four. Below 
we describe two eigenvectors and two generalized eigenvectors corresponding to £ = 0. 

We denote the two eigenvectors lJq, cJq G Ker(£o)- They are given by 



(36) 



where Uq is the NLS solution ([3]). Furthermore, there are two generalized eigenvectors 



9>1 



9,2 



satisfying 



r 1 1^'* 



o;^, ^ = 1,2. (37) 
The two components of the generalized eigenvector ojq'^ (resp. cjq'^) are the first two com- 



ponents of Wg i (resp. Wg 2) given in (31) (resp. (32)). 
4.2 The case e 7^ 

When e is not zero, the algebraic multiplicity of the eigenvalue ^ = is only two and two 
eigenvalues move away from the origin. In order to track them, we use an argument based 
on the Fredholm alternative. 



The linearization of CGL about the T-periodic solution ( 15 ) gives rise to the eigenvalue 
problem (23) defined over the space (L2(T))^ of T-periodic square integrable vector valued 



functions. The linear operator L has two eigenvectors lj^, lj^ G Ker(£). They are given by 

(38) 



'U'' 




' u ' 


U' 




-u 



where U is the solution of CGL given in (15). The eigenvectors (38) can be expanded in e as 



(39) 



where clJq are the two vector- valued functions of Ker(£o) given in (36), and u}\ are given by 

(40) 



where Ui is the coefficient of e in the solution of CGL and is given in (20). Furthermore, the 



vector- valued functions lJq, uj\, and ujI in (39) satisfy the differential equations 



0, 



1,2. 



(41) 
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These equations are obtained by inserting the right-hand-sides of (39) into (23) with i = 
and then writing the relations occurring at the first and second order in e. 

In order to track the eigenvalues of C emerging from the origin, we suppose that lo^ is 
an eigenvector corresponding the a nonzero eigenvalue i^. We expand these in e: 



eh + + 0{e 



(42) 



uj, = uJo + ecji + e^LJz + C(e^). 
The vector is in the kernel of Cq. Thus ojq will be some linear combination of Wq and uJq-. 



0' 



(43) 



where a and P are complex constants. We insert the expressions (42) into the eigenvalue 
problem (23) for C. At first order in e, loi satisfies the following equation: 

Co^^i = h^o - jO-iljq. (44) 



Using (|37|), the first equation in (41) and linear superposition, we obtain the following solu- 

uJi = auj\ + /SljI + ii {aujf + pcjf) , (45) 



tion to (|44[) 



where a and (5 are the constants in (43), u:\ are given in (40), and cjq'* are the generalized 
eigenvectors of Cq satisfying the relations (37). At second order, one finds the following 
equation for UJ2- 

C0OJ2 = (l^J^l + h^^O - jCiLJi - C2^^0- (46) 



We can eliminate the £2 term from (46). Indeed, if we use the second part of (41) and the 

(47) 



expression for uji given in (45), one finds that 



Equation (46) then becomes 



(48) 



We will not solve (48) but rather find a solvability condition that guarantees T-periodic 



solutions. First, we define the usual inner product 



Xi,X2)= / Xnx)X2(x)dx, X, G(L2(r))^ 



(49) 







With respect to (49), the adjoint of Cq from ([t]) is given by 

" d2 



dx^ 



+ 



a 



2Ul 
4f/2 + a 



J, J 



i 
-i 



(50) 



It is then easy to verify that the kernel of £j is generated by J^^q-, i = 1,2, where uJq are the 
generators of the kernel of Cq. 
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We can now obtain compatibility conditions for (48) by taking the scalar product with 
mts of the ker 



elements of the kernel of 



1,2. 



(51) 



Because Jlo^ is in the kernel of Cq, the left-hand side and the last term of the right-hand 
side of (51) are zero. We can then rewrite (51) using the expression for oji given in (45) 



(52) 



+ i2{au:l + (3u;l Ju:^) , z = 1,2. 

This is a homogeneous linear system of equations for a and 13 and we are looking for the 
values of for which (52) has a non-trivial solution. The system can be simplified further 
by making the observation that the components of the vectors Wq, and cJq'^ are even 
functions of x with respect to the axis x = T/2 for i = 1 and odd for i = 2. Thus, any inner 
product in (52) involving two different superscripts will be zero. Furthermore, because J is 
skew-adjoint, the last term of (52) is always zero. We reduce to the two equations 

= ah (£1 {uf, Jul) + {u:l JljI) - (C^ujf, JljI)) , 



(53) 



The system (53) has a non-trivial solution (a,/5) if £1 = 0, which corresponds to the gener 

ith 



ators of the kernel of C given in (38), or if £1 is given by either of the expressions 



or 



(54) 



The components of the vectors uJq and uj\ are real and those of uJq are purely imaginary 
making the expressions in (54) all real. Thus (23) has an eigenvalue on the right side of the 
complex plane for small positive values of e if one of the two solutions in (54) is positive. 
This gives an instability condition for the solution (15). 

Note that in what follows, we make the simplifying assumption that 6 = 1. This is 
done without loss of generality because CGL in the form (13) has the symmetry that if 
q2{x,t) = qi{x,t) = f{x,t) is a solution for a given value of the coefficient r = ro, then 
q2{x,t) = qi{x,t) = 7/(70;, 7^t) is a solution corresponding to the value of the coefficient 
r = rof-y"^. This symmetry can thus be used to set 5 to 1. 



The condition that the first expression in ([54j) be positive gives rise to the inequality 

4fc'2((214fc^ - 223^2 - 15)s + Ak^ + 197A;^ - 225)K^ - 420(1 + s)E^ 
+ 8((214A;6 - 657k^ + 428P + 39)s + Ak^ + 78k^ - 412k^ + 354:)EK^ 

- 168{Ak^ + IQsk'^k'^ + As + 19k'^)E^K^ + 840(s + 2k'^)E^K 

- 210 {l + s)K ((1 + k'^)E^ + (1 + 7k^ + k^)K^ - 2 (1 + 4 /t^ + k^)E K) h 

- 315 k\l + s)K'^ {E -bK k'^K) I2 - 1575 k^ (1 + s) K^h > 0. 

where 



(55) 



Z^(x)dx, I2 



Z^(x) sn^(x)dx, /a 



iK 



Z'^(x) sn^(x)dx. 



12 



(a) 



(b) 



Figure 1: In (a), region 1 (resp. 2) represents the subset of the k-s space in which the 



inequahty (55) (resp. (56)) holds. In (b), region 1 is where the inequality in (16) for the 
existence of the solution of CGL holds. The curve in (b) has an asymptote at A; = and its 
value at /c = 1 is s = — 1/4. In (c), the top curve from (a) and the curve from (b) intersect 
near the point (0.949,0.824); regions 1 and 2 are explained in Figure [2] 



Furthermore, we obtain the following inequality by requiring the second expression to be 
positive 

4A;'2((598fc4 - mSk^ + 799)s - 32k^ - llSfc^ + 169)K^ - 420(1 + s)E^ 
8((598A;'' - 2108fc^ + 2624^;^ - 1138)s - 32k^ + 97k^ + lOAk^ - 193)EK^ 
+ 56((173fc^ - 308A;2 + 143)s + Ik^k'"^ + 8)E'^K^ + 8A0{6sk'^ - 2s + l)E^K 
+ 210 {l + s)K ((1 - 8 - 2 k^)K^ + {Ak^ + 1 k"^ - 2)E K + {I + k'^)E^) h 
+ 315 (1 + s) k^K^ [Sk'^K + E + 3 K) I2 - 1575 (1 + s) k^K^h > 0. 



(56) 



Theorem 3. When either condition (55) or (56) is met, one eigenvalue emerges from the 



origin on the right side of the complex i-space for small positive e. In other words, conditions 



(55) and ^56\ ) are spectral instability criteria for the solution of CGL given in (15). 

Figure [1] shows the regions in the k-s space in which the inequalities for existence (16) 



and instability (55) and (56) are satisfied. Depending on the values of k and s, there are two 



possible types of behaviors for the two emerging eigenvalues (see Figure [2j). In the first case, 
the two eigenvalues move to the left on the real line making our analysis inconclusive since 
other eigenvalues could possibly emerge from zeroes of the Evans function that are elsewhere 
on the imaginary axis. In the second case, the eigenvalues move on opposite directions on 



the real line causing the solution ( 15 ) to be spectrally unstable. 



A Expressions for Vj, V4, Wg, W4, wj, and w 



The first component of Vg is given by 



-6k 



k" dn{6x) sn{5x) Sx + {2k' - 1) ( —bx + Z{bx 



E 
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ReW 



Re(^) 



(a) 



(b) 



Figure 2: Diagram (a) (resp. (b)) represents the behavior of the two emerging eigenvalues 
for values of k and s in region 1 (resp. 2) of Figure [T]^ . 

Furthermore, (Vf)^ = (V^)^ and (Vf).,i = 3,4 are both given by the derivative of (Vg) . 
The first component of V4 is given by 



(v3, = - 



k 

k^ 



E 
K 



k' cn{5x) 5x + cyi{5x)Z{5x) — dn(5a;) sn(5x) 



Furthermore, (Vj)2 = - (V^) ^ and (Vj) . , ^ = 3, 4 are given, respectively, by the derivatives 
of (V4^)^,z = l,2. 

The expressions for Wg, W4, w^, and are rather complicated. However, for the purpose 
of calculating the derivatives of the Evans functions, we only need their expressions evaluated 
at X = T = AK/5: 



-2ikE 
k'^ 



Ak'^E ■ 
6k'^ 

4k^E 
" 6k'^ 

1 

-1 



wl(T) 



2ikE 
Sk'^ 



4k^E 




L k'^ 



w?(T) 



6'^kk'^ 



E{2e - 1) + Kk'^ 
E{2e - 1) + Kk'^ 
26 {k'^K^ + 2 k'^2 e - 1) E K - {3 e - 4k* - 1) 




k'^ 
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2k'^EK + {k^ + l)E^) 

Kk'^ - E 
-Kk'^ + E 

B The Evans Function for Cnoidal NLS Solutions 

In this appendix we use solutions of the AKNS spectral problem for focusing NLS solution 

qo{x,t) = 6ke-''''cn{6x;k), a = 6^1 - 2k^) (57) 

(corresponding to ^) to construct solutions of the ODE system ^ using squared AKNS 
eigenfunctions. We demonstrate that for all but five values of i along the imaginary axis 
(including i = 0), such solutions form a basis for the solution space of ([s]). Using this 
information, we determine all zeros of the Evans function (other than possibly the four other 
imaginary values) using the Floquet spectrum of go- 



-1 



6'^kk'^ 



2 {k'^'K' 



2 {k'^K^ - 2 



B.l Squared Eigenfunctions and the Evans Ansatz 

In this subsection, we will assume that g is a genus one finite-gap solution of focusing NLS 
([T]). These solutions have the form 

g(x,t) = e-^'^*f/o(e), ^ = x-ct, 

for real constants a, c and a T-periodic function Uq which is not necessarily real. (We 
will later specialize to the case of the cnoidal solution above.) The change of variables 
Ui{^,t) = e^"*gi(x,t), U2{C,,t) = e~^"^^q2{x,t) will allow us to linearize about a stationary 
periodic solution when it is applied to system ([T]), yielding 

lUit + - icuif: + aui + 2u\u2 = 0, 

2 

—iU2t + U2^^ + icu2^ + au2 + 2uiU2 = 0. 

(Note that, here, the time derivative holds C, fixed.) We linearize this system about the 
solution ui = = Uo{^) , yielding 

ivit + vi« - icvic + avi + 4|[/o|^Vi + 2UnV2 = 0, 

0—2 (59) 
-iv2t + + icv25 + av2 + 4|t/o| V2 + 2Uo vi = 0. 

We will now show how to obtain solutions for this linearized system, satisfying the ansatz 

vi = eV(0, y2 = e''w2{0- (60) 
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Recall the AKNS system for focusing NLS |2l |3]: 



-iA iq 
iq iA 



i(|gp-2A2) 
2iAg + g^; 



2iAg - 
i(2A2 



l9l 



(61) 



As is well known [121 ET] , the squares of the components of xf) give a solution of the lineariza- 
tion of ([!]) at q. This linearization is 

m + + 4|gp5' + 2q^h = 0, 



-i/it + /;,^^ + 4|g|2/i + 2g2^ = 0, 
and the solution is given hj g = ip\, h = ipl- It is easy to check that the substitutions 



g(x,t) 



-\at 



vi(^,t), h{x,t) 



-.iat 



(62) 



(63) 



turn (l62j) into (jSgj). 

It is also well known that one can use Riemann theta functions to produce formulas 
for finite-gap solutions of NLS and the corresponding solutions of (61), starting with a 
hyperelliptic Riemann surface E of genus g and certain other data (see |3] or [5]). The 
solution (57) arises for g = 1, with the four branch points of E in conjugate pairs in the 
complex plane, so that the equation of S is 



(A-Ai)(A-Ai)(A-A2)(A-A2), 



(64) 



where ReAi < ReA2. (See §3 of [3] for a derivation of these solutions from general finite- 
gap formulas.) In genus one, the finite-gap solutions of (61), which are known as Baker 
eigenfunctions, take the form 

iji = exp (i(fii(P) - f )x + i(fi2(P) + f )t) ei(e, P), 

V^2 = exp {m^{p) + f )x + i(^]2(P) - f )t) e2(e, P), 

where E,N are constants determined by S, P = (A,/i) is an arbitrary point on S that 
projects to A e C, ^1^2 are certain Abelian integrals on E, and 6i_2 will be defined below. 
For the moment, we note that 0i, 02 depend only on P and = x — ct (where c is determined 
by the branch points of S) and are T-periodic in ^. 

After the substitutions (63) the squared Baker eigenfunctions yield solutions to (59) of 
the form 









72. 




e-'"V2 



,2i5ni (P)+2it{n2(P)+cQ.i (P)) 



(^{cE-N-a)t+\EiiQ^(^^^ p)2 



In particular, these vi, V2 satisfy the ansatz (|60j) for a = cE — N and 

£ = 2i(fi2(P) + cl]i(P)), 

with 

w,{0 = e'(2^i(^)-^)«0i(e; P)^ W2iO = e^(2^^(^)+^)«e2(e; P)'- 
Proposition 4. The formula (65) implies that i = 4i/i. 



(65) 
(66) 
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Proof. The surface S has genus 1, and a basis a and b (shown in Figure [3]) for homology 
cycles. The differentials df2i,df22 have zero a-period. Denoting their 6-periods as V,W 
respectively, from [3] we have 

W/V = -c = (Ai + A2 + AT + A^). 

(Note that the meaning of c here is —1 times its definition in ^.) Thus, the differential 

d£ = 2i(dfi2 + cdfii) 

has zero a- and 6-periods. So, although ^1^2 sue not well-defined on S, and so £ is a well- 
defined function on S. 

The integrals fli use the branchpoint A2 as basepoint, so ^ = there. Because L*dfli = 
—dfli for i = 1,2 (where l is the sheet exchange automorphism (A,/i) t— > (A, — yu)), then 
L*i = —i. So, it follows that i vanishes at all four branch points, which are the same points 
at which the coordinate /i vanishes. The surface S has two points at infinity, called oo±, 
where both A and /x approach complex infinity and 

A^^-l^-^(^n^" + ''(^"') 
respectively. Combining the asymptotic expansions for fli and ^2 gives 

^ = ±2i (2 + cA-^ + C»(A-^)). 

Therefore, 

= 4i + C(A-2). 

Because i/fi is bounded and holomorphic on S, then it must be equal to the constant 4i, 
i.e., 

i = 4i/i. 

In other words, we have just proved the identity 

n2{P) = -cni{P) + 2fi (67) 
for genus one NLS solutions. □ 



B.2 Forming a Basis of Solutions 



Proposition 4 and (64) imply that for generic values of i there are four distinct points P G S 



for which the Baker eigenfunctions may be used to construct a solution of the eigenvalue 
problem The exceptional values are those for which two roots of (64), as a polynomial 
equation for A, coincide. 

From here on, we specialize to the case where the branch points in (64) satisfy A2 = — Ai, 
with ReAi 7^ 0, which yields the NLS solution (57). In this case, c = and E = V/2, the 

—k' + ik, and the finite-gap 

2|Ai| 



branch points are related to the elliptic modulus by Ai/|Ai| : 
solution coincides with (57) with 6 = 2|Ai|. The exceptional values of /i are 



fj, = ±|Ai| 



fi = ± Im(A^), 



(68) 
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corresponding respectively to a double root occurring at A = 0, or two pairs of roots coin- 
ciding at opposite points along the real axis (if k"^ < 1/2) or imaginary axis (if k"^ > 1/2). 
When /i is not an exceptional value, we will construct a matrix solution W{^;i) of J8) 



where each column is of the form [wi{^),W2{^),w[{^),W2{^)] and Wi,W2 are given by (66) 



for one of the four points P = (A, /x) satisfying (64) for the given /x. For purposes of showing 
that the matrix solution is nonsingular, it suffices to evaluate its determinant at one value 

ofe 



We now need to specify the form of the factors 0i, G2 in (66). Specializing the formulas 
given in §5 of [5], and using ^ = x — ct = x, we have 

^ eiAiP) + iVx-D)eiD) ^ ,p, e{A{p) + iVx-D-R)e{D) 

e{iVx-D)e{A{p)-Dy ^^i^,^; e{iVx-D)e{A{p)-D) ' 

In these formulas, 

m V = 7r|Ai - \2\/K = 27r\Xi\/K; 

• A{P) is the Abel map, obtained by integrating the differential oo = VdX/{2fi) from 
oo_ to P on E; 

• ^^3(P) is the integral, from A2 to P, of the unique meromorphic differential dfl^ on P 
which has zero a-period and satisfies dAi ~ ±A~^dA near oo±; 



R is minus the period of d^a on the cycle b running from Ai to A2 (see Figure [3]), and 
is also the value ^4(00+), when the path of integration is chosen to avoid the homology 
cycles a and b (see §5 in [_5J); in the case when c = 0, R = tt{K' /K + i) (see §3 in |5J); 

D is an arbitrary constant which is pure imaginary; 

6{z) is the Riemann theta function with period 27ri and quasiperiod B = —2'kK'/K. 
It is related to the Jacobi theta functions of modulus k by 6{z) = 93{z/ {2i)). 




Figure 3: Homology cycles and other integration paths on the genus one Riemann surface 
with complex conjugate branch points. Solid curves appear on the upper sheet, dashed curves 
on the lower sheet, and branch cuts extend between each branch point and its conjugate. 
The horizontal line represents a set of points where n is real and positive. 
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In genus one, the constant D has no significance, as it can be absorbed through a shift 
in X. Thus, we may assume D = in the above formulas. Once this is done, the columns of 
for ^ = take the form 



(69) 



'wi{0)' 




1 


W2{0) 






^(o) 




2ig{P) 
-2if{Pfh{P) 


^(0) 





where 



e{A{p) 



0{A{P)) 
h{P) := Qi{P) + VQ{A{P) 



g{P) ■.= n,{P) + VQ{A{P))- 



E 



E 
2"' 



and Q{z) = 6'{z)/6{z) is the logarithmic derivative of the Riemann theta function. (Note 
that A{P) and fli{P) are not individually well-defined on E, because of the nonzero periods 
of the corresponding differentials along the homology cycles. We take the convention that 
the paths of integration for these differ by a fixed path T (shown in Figure |3| from oo_ to 
A2 in Eq, where Sq denotes the simply-connected domain that results from cutting S along 
the homology cycles. With this convention, /, g, h are well-defined meromorphic functions 
on E.) 

The Riemann surface E has two holomorphic involutions, namely l : (A,/i) 1— >■ (A, — /i) 
and a : (A,/i) (— A,/x). Because L*dfli = —dQi, l*uj = —u and l fixes the basepoint for 
i = 1, 2, 3, it follows that 



Because a*dQi 



f{iP) = 1//(P), giiP) = -h{P). 
-dfii, a*dfi3 = dfis — u;, a*uj = —lo, and E = V/2, it follows that 
fiaP) = /(P), giaP) = -g{P). 



Let Pi,P2,P3, Pi be four points on E corresponding to a given (non-exceptional) value of 
/i. Without loss of generality, we may assume that P2 = cr(Pi) and P4 = criPs). Using the 
above formulas for the behaviour of /, g, h under a, we find that the matrix with columns 
given by (69) is 



W{0) 



-f{Pif 
2i^7(Pi) 
-2\f{P,fh{P,] 



1 

-f\Pi? 
-2i(7(Pi) 
2if{P,fh{P,) 



1 

-f{P^f 



1 

-f\P^? 



2i^(P3) -2ig{P3) 
-2i/(P3)2/i(P3) 2if{P,Yh{Ps) 



Proposition 5. The functions g, h satisfy the identity 

-g{P) = h{P) = X, P=(A,/i)GE. 

Proof. The poles of g{P) = ^i{P) + Q{A{P)) — E/2 can only occur where Qi has a pole 
or 9{A{P)) = 0. (We now work on the cut Riemann surface Eg.) These occur only at oo_ 
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(near which ^]l ~ -A + C(l), but 6{A{P)) = 6(0) = ^IKj-n ^ 0), and at 00+ (near which 
f2i ~ A + 0{\) and 6'(A(P)) = Q[K) = 0). Near oo_,_, use w = as local coordinate. Then 
A{P) = ^ /i~^dA implies that 

dA{P) V V 2, 



and therefore A{P) = R — + 0{w^). Because 9{z) has a simple zero at 2; = i?, 

giP) = n,iP) + ^^^^ -^ = w-' + V + Oil) = -A + 0(1). 

The oddness of g with respect to a implies that g{P) = when P lies over the origin in 
the A-plane. Therefore, the quotient g{P)/X is a bounded holomorphic function on S, and 
so must be a constant. The above asymptotic expansion for g implies that g{P) = —A, and 
it follows from the fact that /i(tP) = -fi'(P) that h{P) = A. □ 

Taking Prop. |5]into account, we compute 

detiy(O) = -16A(Pi)A(P3)(/(Pi)' - fiPsry. (70) 



Proposition 6. For /i not equal to zero or any of the exceptional values in (68), the matrix N 
is nonsingular. (We conjecture that the analogous statement is true for genus one finite-gap 
solutions in general.) 



Proof. Because we are excluding the exceptional values, neither of the A- values in (70) are 
zero, and it only remains to establish that /(Pi) 7^ ±f{Ps). 

Define F{P) = f{P) + 1//(P), which is a well-defined meromorphic function of A. The 
formula for /(P) shows that its only pole is a second-order pole at 00+ (because exp(f23(P)) ~ 
A near there), and hence P is a second-order polynomial in A. It is known that /(±Ai) = 1 
and /(±Ai) = —1 (see §7 of [6]). Hence, 

P(A) = i^^^, d:=Re{Xl). 
Ai — Ai 

Using the quadratic formula, and the fact that / has a second-order zero at oo_, we obtain 

AP) - '-^If^. (71) 
Ai — Ai 

Let yi = A(Pi) and y^ = XiP-s)- Because yi and y^ are distinct roots of 

{X'-Xl){X'-T^')=fi', (72) 

and yi 7^ —y^, then yi + yi = 2d. Hence, P(yi) = 2{yf — yi) = —F^y^), and therefore 
/(Pi) 7^ /(-Ps)- Furthermore, using (71), we see that 

(A? - AT') (/(Pi) + /(P3)) = 2{yl + yi -2d + 2/i) = 4/.. 
Thus, because we assume 7^ 0, /(Pi) 7^ — /(P3). □ 
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Assuming that /x is not zero or one of the exceptional values, we can define the transfer 
matrix 

N = W{T)W{0)-\ 

whose eigenvalues describe the growth of solutions to (|8| over one period. (In particular, 
there is a periodic solution if and only if has eigenvalue one.) Using the fact that T = 
An/V = 2tt/ E and VLi{a P) = —Qi{P) + E/2 on the cut surface, we calculate that 



W{T) = W{0) 



g2ini(Pi)T Q Q 

Q g-2iQl(Pl)T Q 

Q Q g2iQl(P3)T 









g-2if7l(P3)T 



Then the characteristic polynomial of the transfer matrix is 

det{N~Tl) = ((r+ 1)2 -4rcos(^]l(P)^))((r + 1)2 -4rcos(fii(P3)r)). 
By substituting r = 1 in this formula, we obtain the 



Proposition 7. For fi not equal to zero or the exceptional values in (68), the system (|8]), 
where i = 4i/i, has a periodic solution if and only if sm{fli{Pi)T) = or sin(f2i(P3)T) = 0. 



B.3 Imaginary Zeros of the Evans Function 

For periodic finite-gap solutions of NLS, the Floquet discriminant [5j takes the form 



'2tt 

A(A) = 2cos ( yn^iP) 



2cos( |l^i(P) 



(While Qi is odd with respect to the involution l, the even-ness of the cosine makes A a 
well-defined function of A.) This has the property that the AKNS system (61) admits a 
T-periodic (respectively, antiperiodic) solution if and only A = ±2. 

The consequence of Propositions |6] and [7] is that, for values of i that correspond to 



nonzero values of excluding those in (68), the Evans function is equal to zero if and only 
if A = ±2 or A = at an opposite pair out of the four A-values corresponding to /x. When 
this happens, the corresponding pair of columns of W are periodic, and the corresponding 
value of £ is a zero of the Evans function of geometric multiplicity two. Thus, we can use 
the Floquet discriminant to find the (nonexceptional) zeros of the Evans functionj^ 
In §3.5 of 0, the discriminant for genus one finite-gap solutions is calculated as 



A(A) = -2 cos 



2iK Z(u) - p' 



, cn M dn M sn u 
1 — /S^ sn^ u 



where the variable u is related to A by 



sn u 



V.(A) := 



(Ai-A2)(A-Ai) 
(Ai-Ai)(A-A2) 



(73) 



(74) 



^Even if one of the double points gives an exceptional value of /i, it still gives a zero of the Evans 
function. However, there may be additional zeros at exceptional values of fj,, corresponding to possible 



periodic solutions of i&h obtained from solutions of the form ( 66 1 by reduction of order 



21 



and the parameter (3 and the modulus k are related to the branch points by 



2 ^(Ai — Ai)(A2 — A2) „2 ^1 ~ -^1 



= — -^^^ — ^, 

|Ai — A2P Ai — A2 

In the symmetric case (i.e., A2 = — Ai) which we are considering here, these specialize to 

-A;' + iA; = Ai/IAil, = k{k - ik'), 

and the variable u is determined by A as follows. The function sn^ u has periods 2K and 2iK', 
and is even in u, so it suffices to restrict u to the set (0, K] x {—K', K'] in the complex plane. 
Then there is a unique u in this set for each A in the extended complex plane. Conversely, 



solving (74) for A in terms of u gives 

Ai sn^ M + Ai/?~2 



X{u) 



^ 2 _ g^2 y 



The locus of A- values for which A (A) is real and between 2 and —2 is known as the 
continuous Floquet spectrum of the NLS potential q. For genus one solutions, it consists of 
the real axis and two bands terminating at the branch points. (For low values of k the bands 
emerge from the real axis, but for k sufficiently near 1, the bands become detached from the 
real axis (see Figure |4]); the transition occurs around k ~ .9089 |14j.) 

The periodic points (i.e., where A = ±2) consist of the branch points and countably 
many points (known as double points) on the real axis. Points where A = will naturally 
interlace the periodic points along the real axis, but will also occur along the imaginary axis 
if the bands are detached. When this happens, though, the corresponding value of /x is real. 
Hence, all the zeros of the Evans function will occur along the imaginary axis in the complex 
£-plane. 

Figures |4] and [5] show the level sets of A(A) for several different moduli, and the zeros of 
Evans function in the £-plane. The latter figures also include the locus of £- values which are 
related by ([72|, with i = 4i/i, to A- values in the continuous spectrum. 
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Figure 4: Level curves of the Floquet discriminant A in the complex A-plane for elliptic 
moduli k = .6, k = .8 and k = .92, respectively. (We take |Ai| = 6/2 = 1.) In each of the 
diagrams, the black curves indicate where A is real, the gray curves where A is imaginary, 
the dashed curves where Re A = ±2, and the dotted curves where is real (which includes 
the real and imaginary axes). Thus, branch points occur where dotted intersects dashed, 
periodic points occur where black intersects dashed, and points where A = occur where 
black intersects grey. 
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Figure 5: Zeros of the Evans function in the complex £-plane, for k — .6, k — .8 and k — .92, 
respectively. Note that all zeros are of multiplicity two, except for the origin, which is of 
multiplicity four. The continuous spectrum consists of the imaginary axis and the figure 
eight curve. 
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